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Abstract. Rotation appears as a dominant effect in massive star evolution. It 
largely affects all the model outputs: inner structure, tracks, lifetimes, isochrones, 
surface compositions, blue to red supergiant ratios, etc. At lower metallicities, 
the effects of rotational mixing are larger; also, more stars may reach critical 
velocity, even if the initial distribution of rotational velocities is the same. 



1. Introduction 

In this review, we shall concentrate mostly on the rotational properties of mas- 
sive stars. Rotation has an importance which for massive stars in the Galaxy 
is comparable to that of mass loss by stellar winds. For massive stars at lower 
metallicities, it is likely that rotation plays an even more important role. 

2. Structure, mixing and mass loss of rotating stars 

2.1. Structure 

The equations of stellar structure for a rotating star were written by Kippenhahn 
& Thomas (1970). These equations apply to a star in solid body or cylindrical 
rotation. For stars in differential rotation, the use of the above equations is 
incorrect. Here, we consider the case of shellular rotation with $7 = J7(r) (Zahn 
1992). In this case, the structure equations need to be modified (Meynet & 
Maeder 1997). Structural effects due to the centrifugal force are in general small 
in the interior. However, the distorsion of stellar surface may be large enough 
to produce significant shifts in the HR diagram (Maeder & Peytremann 1970). 

2.2. Internal transport 

For the transport of chemical elements, we consider the effects of shear mixing, 
of meridional circulation and their interactions with the horizontal turbulence. 
For shellular rotation, the equation of transport of angular momentum in the 
vertical direction is in lagrangian coordinates (Zahn 1992), 
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J7(r) is the mean angular velocity at level r. U{r) is the radial term of the 
vertical component of the velocity of the meridional circulation. £> s hear is the 
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coefficient of shear diffusion. From the two terms in the second member, we see 
that advection and diffusion are not the same. Thus, the transport of angular 
momentum by circulation cannot be treated as a diffusion. For the changes 
of the chemical elements due to transport, we may use a diffusion equation 
with a diffusion coefficient which is the sum of .D s hear> as given below, and of 

30Dh ■ Dc{[ ex P resses tne resulting effect of meridional circulation and 
of a large horizontal turbulence (Chaboyer & Zahn 1992). This expression of 
D e Q tells us that the vertical advection of chemical elements is inhibited by 
the strong horizontal turbulence characterized by D^. The usual estimate of 
Z?h = |2V(r) — aU(r)\ was given by Zahn (1992). A recent study suggests 
that this coefficient is at least an order of magnitude larger (Maeder 2002), 
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For n=l, 3 or 5 A « 0.134, 0.0927, 0.0782 respectively. For the coefficient D shcar , 
we use the results by Maeder (1997), Talon & Zahn (1997) and Maeder k Meynet 
(2001). -Dshcar is also modified by the horizontal turbulence Z?h, in the sense that 
a larger leads to a decrease of the mixing of elements, 



-^sncar — r 1 r 



|v^(i + ^) + (v ad -v rad ) 



«(>«iE7)'-^-*> 



with b=0.8836. The radial component U(r) in Eqs. (1) and (2) was given by 
Maeder & Zahn (1998), 
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where = M (l — 2 ^G Pm ) i s * ne reduced mass, with the notations given in the 
quoted paper. The driving term in the square brackets in the second member is 
Eq. It behaves mainly like, 
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The term 7> means the average on the considered equipotential. The term with 
the minus sign in the square bracket is the Gratton-Opik term, which becomes 
important in the outer layers due to the decrease of the local density. It can 
produce negative values of U(r). A negative U (r) means a circulation going down 
along the polar axis and up in the equatorial plane. This makes an outwards 
transport of angular momentum, while a positive U(r) gives an inward transport 
of angular momentum. Before, the term V M was never properly included. The 
justification for this term is not so straightforward (Maeder & Zahn 1998). 

2.3. Critical velocity and mass loss in a rotating star 

For a rotating star, one must consider the flux F(i!}) at a given colatitude $ as 
given by von Zeipel's theorem. This theorem has been studied for a star with 
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Figure 1. Left: The second critical velocity v cr it.2 expressed as a fraction of 
«crit,i plotted as a function of the Eddington factor r max , which is the largest 
value of the Eddington factor over the stellar surface. Right: Evolution of the 
surface velocities for a 60 M Q star with 3 different initial metallicities. 



shellular rotation by Maeder (1999), who finds F($) = - ^gmJ ckII + ((#)], 
which shows a small extra-term depending on CI, T, VT, p and opacities. 
In a rotating star, the Eddington factor becomes a local quantity T^(i9). We 
define it as the ratio of the local flux F($) given by the von Zeipel theorem to 
the maximum possible local flux, which is Fy im ($) = — ^ySeff ($)• Thus, one has 

where the opacity depends on the colatitude since T e s also depends 
The Eddington factor depends on the angular velocity O on the isobaric surface. 
This shows that the maximum luminosity of a rotating star is decreased by 
rotation. It is to be stressed that if the limit To^i?) = 1 happens to be met in 
general at the equator, it is not because g c g is the lowest there, but because 
the opacity is the highest! Indeed, the dependences of F({)) and of -Fii m ($) with 
respect to g e g have cancelled each other in rn(#). 

Often, the critical velocity in a rotating star is written as v^ lit = ^-(1 — T). 
This expression is incorrect, since it would apply only to uniformly bright stars. 
The critical velocity of a rotating star is given by the zero of the equation 
expressing the total gravity gt ot = fiCfr + ftad = Sgmv + 5rot + ftad- This is 

Stot = 02r[i-r n (0)] • (7) 

This equation has two roots (Maeder & Meynet 2000). The first that is met 

determines the critical velocity. The first root is as usual f cr it,i = (|§^) 2 > 
where i? p b is the polar radius at break-up. The second root applies to Eddington 
factors bigger than 0.639. It is illustrated in Fig. 1 (left), which shows that, 
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Figure 2. Left: Iso-mass loss distribution for a 120 M Q star with 
log -jf^ = 6.0 and T g = 30000 K rotating at a fraction 0.80 of break-up ve- 
locity. Right: The same with T cff = 25000 K, (Maeder & Desjacques 2000). 



when r is close to 1.0, the critical velocity goes to zero. The theory of radiative 
winds applied to a rotating star leads to an expression (Maeder 1999) of the 
mass flux as a function of colatitude. Figs. 2 illustrate the distribution of the 
mass loss rates around a massive star of 120 M with two different T e g. For 
a star hot enough to have electron scattering opacity as the dominant opacity 
source from pole to equator, the iso-mass loss curve has a peanut-like shape 
(Fig. 2 left). This results from the fact that the pole is hotter ( "<7 e fj-effect" ) . If 
the T e ff of the star is lower, a bistability limit (i.e. a steep increase of the opacity, 
cf. Lamers 1995) may occur somewhere between the pole and the equator, due 
to the decrase of T c g from pole to equator. This "opacity-effect" produces 
an equatorial enhancement of the mass loss (Fig. 2 right). The anisotropies 
of mass loss influence the loss of angular momentum, in particular polar mass 
loss removes mass but relatively little angular momentum. This may strongly 
influence the evolution (Maeder 2002; see Fig. 5 left below). We may estimate 
the mass loss rates of a rotating star compared to that of a non-rotating star at 
the same location in the HR diagram. The result is (Maeder & Meynet 2000) 

M(o) ^ a-r)^ 1 (8) 
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where T is the electron scattering opacity for a non-rotating star with the same 
mass and luminosity, a is a force multiplier (Lamers et al. 1995). Values of 
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Figure 3. Left: Evolution during the MS phase of the angular velocity in 
the interior of a 20 M Q star at Z=0.02. Right: the same for a 15 M Q star 
at Z = 10~ 5 . The values of the central H-content X c are indicated. (From 
Meynet & Maeder 2000, 2002). 



a are given in Table 1, which shows, for stars of different masses and T, the 
values of the ratio of the mass loss rates for a star at critical rotation to that 
of a non-rotating star of the same M, L and T e g. The values oo indicate that 
the critical velocity t> cr it,2 is met. Close to F = 1 the increase of the mass loss 
rates may be quite large. For blue and red supergiants, the increase could also 
be large, however supergiants do not rotate very fast in general. 

We must distinguish 3 cases of stellar break-up: 1.— The T-Limit, when 
radiation effects largely dominate; 2.- The fi-Limit, when rotation effects are 
essentially determining break-up and 3.- The fiT-Limit, when both rotation 
and radiation are important for the critical velocity and Eq. (8) applies. For 
very massive stars, the Or-Limit is clearly the relevant case. 



Table 1. Values of ■ , ( for different stellar masses. The four force multi- 

M(0) 

pliers a apply respectively to stars with: 4.70 > logT c ff > 4.35 (type B1.5 or 
earlier), logT eff = 4.30 (type B2.5), 4.00 (B9.5), 3.90 (A7) respectively. 
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Figure 4. Evolution of the equatorial velocities along the redward parts of 
the tracks in the HR diagram for initial velocities of 300 km s _1 and Z = 0.02, 
(Meynet & Maeder 2000). 



3. Evolution of rotational velocities 

Fig. 3 shows the evolution during the MS phase of the internal angular velocities 
in models at Z = 0.02 (left) and at Z = 10 -5 (right). We see that the decrease 
of the core rotation is stronger at higher Z. This is due to the larger mass 
loss which removes more angular momentum, (the case of isotropic mass loss 
is considered here with Eq. (8) above). We note that the gradient of Q built 
during MS evolution outside the core is larger at lower Z. There are 2 reasons 
for that. One is the higher compactness of the star at lower Z. The second one is 
more subtle. At Z = 0.02, the density of the outer layers is lower and the star is 
more extended than at lower Z, thus the Gratton-Opik term is more important. 
This produces an outward transport of angular momentum and smoothens the 
O-gradient at Z = 0.02 more than at lower Z. The steeper ^-gradient at lower 
Z is leading to stronger mixing of chemical elements, as shown below. 

The evolution of the surface velocities results from internal transport and 
surface losses of angular momentum. Fig. 1 (right) shows the evolution of the 
surface velocities for a 60 M model for 3 very different metallicities Z. For such 
a large mass, the high losses of mass and angular momentum lead to a strong 
decrease of rotation at Z = 0.02. On the contrary, at lower Z the velocities are 
growing, as a result of the contraction of the stellar core and of the relatively 
good coupling of the outer layers produced by circulation. Thus, it is likely 
that at low Z, rotation may be a dominant effect in massive star evolution. 
These stars may reach break-up and then loose a significant amount of matter. 
Paradoxically, low Z stars may loose mass due to fast rotation as a consequence 
of their initial slow winds! 
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Figure 5. Left: Evolution of the ratio tj- of the angular velocity to the 
critical value during the MS phase of a 40 M Q star with an initial velocity of 
500 km s" 1 (Maeder 2002). Right: Evolutionary tracks on the MS for 20 M 
models of different Z and initial velocities v- m i, (Meynet & Maeder 2002). 



In addition, we may wonder whether the initial distribution of stellar rota- 
tion velocities is the same at lower Z. There are some reasons (Maeder et al. 
1999) to suspect that the initial rotation velocities may be faster at lower Z. 

Fig. 4 shows the lines of iso-rotational velocities in the HR diagram of 
massive stars at Z = 0.02 for an initial velocity of 300 km/s. We see an important 
decrease for the most massive stars. This provides the basis for interesting 
comparison tests. Let us emphasize that the decline will be less important when 
the anisotropic losses of angular momentum are accounted for, which is not the 
case in Fig. 4. Fig. 5 (left) shows the very different evolution of the surface 
rotation of fast rotating stars for different cases of mass loss. For zero-mass 
loss, critical velocities are reached very quickly. For isotropic mass loss (here 
the relatively lower rates by Vink et al. 2001), there is no increase of j?-. The 
case with account of the anisotropic losses of angular momentum gives results 
intermediate between the two previous cases, with some significant increase of 
jy-. If there is some magnetic coupling, jj- will be larger during evolution. Thus, 
the study of rotational velocities for stars at various distances from the ZAMS 
may provide a powerful test on the role of magnetic field. 



4. HR diagram, lifetimes and isochrones 

Table 2 provides various information on the effects of rotation for a 20 M 
star. We show the average velocities during the MS phase, the MS lifetimes, the 
masses and velocities at the end of the MS, the helium content, the N/C and 
N/O ratios at the end of the MS. We notice the growth of the MS lifetimes for 
higher rotation velocities, as well as the lower final masses. This results from 
both the longer MS lifetimes and from the higher mass loss rates. We notice 
the chemical enrichments at the stellar surface. Some isochrones calculated for 
initial velocities Vi n \ = 200 km/s are displayed in Fig. 6 (left); these Vi n i lead to 
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Figure 6. Left: Isochrones with and without rotation for solar metallicity. 
Right: Evolutionary tracks in the log g c ff vs. log T c ff diagram. The dashed 
lines correspond to non-rotating tracks, the continuous lines are for models 
with an initial rotation of 300 km s _1 . The small dotted line corresponds to a 
very fast rotating 60 M©. The observations are by Herrero et al. (2000) and 
previous works. 

average velocities on the MS of about 140 km/s. If we assign cluster ages from 
these isochrones, we obtain ages typically 25% larger than from the standard 
models without rotation. For average velocities of about 220 km/s, the difference 
on the age estimate would be larger. This effect may easily help us to reconcile 
the ages for the Pleiades determined from the turnoff with the ages from lithium 
abundances in the low mass stars of the Pleiades (Martin et al. 1998). 

Table 2. Properties of 20 M Q models at the end of the MS for different 
initial velocities. The velocities are in km s _1 , the lifetimes in million years, 
the masses in solar mass and the surface abundances in mass fraction. 



Vini 


V 


in 


Al 


V 


He 


N/C 


N/O 








7.350 


19.019 





0.30 


0.25 


0.12 


50 


30 


7.720 


18.896 


18 


0.30 


0.27 


0.12 


100 


62 


8.292 


18.681 


16 


0.30 


0.45 


0.19 


200 


132 


8.901 


18.324 


91 


0.32 


1.01 


0.38 


300 


197 


9.309 


18.020 


167 


0.35 


1.77 


0.58 


400 


253 


9.745 


17.646 


217 


0.37 


2.54 


0.76 


500 


294 


10.275 


17.181 


213 


0.40 


3.65 


0.99 


580 


304 


10.324 


17.148 


214 


0.39 


3.75 


1.00 



Fig. 6 (right) shows the tracks in the logoff vs. logT e g- plot.We see that if 
we assign an evolutionary mass to a rotating star on the basis of a track without 
rotation, we obtain a mass which may be up to a factor of 2 too large. Appropri- 
ate tracks have to be used for rotating stars. The non-respect of this prescription 
is likely the explanation of the so-called problem of the mass discrepancy. 
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Figure 7. Left: Evolution as a function of logT c s of the abundance ratio N/C 
(in number), normalized to their initial values. The tracks are for Z = 0.004. 
A model at Z=0.02 is also shown. The track with km/s applies to both 
Z = 0.020 and Z = 0.004. Right: Evolution of T eS as a function of the 
fraction of lifetime spent in the He-burning phase fo 15, 20 and 25 M Q at 
Z = 0.004 for v mi = and 300 km/s, (Maeder & Meynet 2001). 



5. Surface abundances 

As a result of internal mixing, He and N are brought to the surface, while C 
is depleted. At solar Z, the N/C ratios are increased by a factor 2-3 for a 20 
M Q model (Fig. 7 left). For lower Z, we see that the enrichment in N/C is 
much larger. This results from the larger shears in low metallicity models (Fig. 
3). The higher N-enhancements at lower Z are consistent with the results by 
Venn (1998, 1999), who has found excesses of N/H up to about a factor of 10 
in supergiants of the SMC. At Z = 0.004, the models show no sign of primary 
nitrogen. Models at Z = 10 -5 by Meynet <fe Maeder (2002) show large amounts 
of primary nitrogen, which change the sum of CNO elements in advanced stages. 

There is still a problem (Herrero et al. 2000). There are many fast rotating 
stars, which have He excesses. The models indicate that when the products 
of CNO burning appear at the surface, the velocities have already significantly 
decreased. A possible explanation is that there are initially many very fast 
rotating stars, where the enrichments appear quickly. Also, nuclear reactions 
in massive stars are active in the pre-MS phase and thus mixing may already 
occur during the pre-MS phase. 



6. The blue to red supergiant ratio 

In clusters of low Z, as in the SMC, there are lots of red supergiants. However, 
the models of massive stars at low Z do not reach in general the red stage. 
If they do it, this happens only very late in evolution, so that there are no red 
supergiants predicted. Fig. 7 (right) illustrates the point: for zero rotation most 
of the He-burning phase is spent in the blue and very little in the red. 
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Rotation very much changes the above results. Fig. 7 (right) shows that 
the models of 15 to 25 M with rotation spend at least half of the He-burning 
lifetime in the red. This is in agreement with the blue to red supergiant ratios of 
0.5-0.8 observed in the SMC. The physical reason why rotation makes more red 
supergiants is the following one: with rotation, the He-burning core is larger 
and there are also more He mixed just outside the core, in particular at the 
location of the H-burning shell. Since the core is larger, and as there is also 
more He outside, the H-shell is less active and the opacity is also lower. Thus, 
while there is an intermediate convective zone in non-rotating models, there is 
none in models with rotation. Convection means a polytropic index n = 1.5 and 
thus a high degree of compactness, which keeps the non-rotating stars in the 
blue. With rotation, the absence of intermediate convective zone permits the 
envelope to largely expand and the star moves toward the red supergiant stage. 

On the whole, we see that rotation is an essential ingredient of stellar mod- 
els and that rotation modifies all the model outputs. 

Acknowledgements: We express our thanks to Raphael Hirschi for a careful 
reading of the manuscript. 
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